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Abstract 

We use the AdS/CFT correspondence to show that the heavy quark (static) potential in a 
strongly-coupled plasma develops an imaginary part at finite temperature. Thus, deeply bound 
heavy quarkonia states acquire a small nonzero thermal width when the t'Hooft coupling A = 
g 2 N c 3> 1 and the number of colors N c —* oo. In the dual gravity description, this imaginary 
contribution comes from thermal fluctuations around the bottom of the classical sagging string in 
the bulk that connects the heavy quarks located at the boundary. We predict a strong suppression 
of T's in heavy-ion collisions and discuss how this may be used to estimate the initial temperature. 

PACS numbers: 12.38.Mh, ll.25.Tq, ll.25.Wx, 24.85. +p, 12.39.Pn 



1 



The conjectured equivalence of strongly-coupled 4-dimensiona. N = 4 Supersymmetric 
Yang-Mills (SYM) to type IIB string theory on AdSs ® S§ [1] has led to new insight into 
the strong coupling dynamics of large N c gauge theories at finite temperature. In fact, the 
Anti-de Sitter/Conformal Field Theory (AdS/CFT) correspondence has been particularly 
useful to compute real time correlators of gauge invariant quantities in strongly-coupled 
olasmas such as 2-point functions of the energy-momentum tensor at finite temperature 
2J. For instance, it was shown that the shear viscosity to entropy density ratio satisfies 
in rj/s > 1 in all strongly-coupled gauge theories that possess a dual description in terms 

n 

of supergravity [3[. Here, we adopt strongly coupled M = 4 SYM as a toy model for the 
deconfined, high temperature phase of QCD. 

The dual description of the gauge theory at finite temperature involves a near-extremal 
black brane in the bulk, which leads to a 5-dimensional metric (in real time) given by 

ds 2 = -G m (U)dt 2 + G xx (U)dx 2 + Guu{U)dU 2 . (1) 

Goo(Uh) = defines the location of the black brane horizon in the 5th coordinate, and 
the boundary is at U — > oo. 

The potential between fundamental static sources separated by a distance L at large t' 
Hooft coupling A in J\f = 4 SYM was computed in {4, 5] and shown to be proportional to 
1/L (due to conformal invariance of the theory) and to vA, which indicates that charges are 



partially screened even in vacuum 



HE 



6|. In general, thermal effects are expected to reduce 
the binding energies of small states of very heavy quarks at high T. At strong coupling 
thermal screening corrections appear at the same order in A as the vacuum potential 
as opposed to Debye screening in weakly coupled quark-gluon plasmas [8j. However, at 
distances L < 1/T these corrections are suppressed by a factor of (LT) A , which originates 
from the behavior of the dual bulk geometry near the black brane horizon. Thermal effects 
also diminish as rj/s increases jo| . 

In this Letter we show that at finite temperature the static potential in a strongly-coupled 
plasma develops an imaginary part due to fluctuations about the extremal configuration, 
which corresponds to a string connecting the fundamental sources at the boundary of the 
geometry. Such an imaginary part arises also in perturbative Quantum Chromodynamics 
(pQCD) at order ~ g 4 due to Landau damping of the static gluon exchanged by the heavy 
quark sources |10 |. At large t' Hooft coupling, however, it appears already at the same order 
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as the vacuum potential, i.e., at 0(vA). Therefore, energy levels in this potential are not 
sharp because they acquire a thermal width, E = E v3iC + AEt — iT. The width T is smaller 
than the vacuum energy E vac if LT < 1 (which is the relevant regime in the limit of very 
heavy quarks, niq — > oo) and of the same order in both A and LT as the shift AEt of the 
real part of the potential due to thermal screening effects. We propose that the imaginary 
part of the potential mentioned above can be observed experimentally via the suppression 
of T to dilepton decays in heavy-ion versus p+p collisions at RHIC and LHC. 

The relevant operator for our discussion is the path-ordered Wilson loop defined as 

W(C) = J-TrPe^^" (2) 

where C denotes a closed loop in the boundary, is the non-Abelian gauge field, and the 
trace is over the fundamental representation of SU(N C ). We consider a rectangular loop with 
one direction along the time coordinate t and spatial extension L. In the asymptotic limit 
t — > oo, the vacuum expectation value of the loop defines a static potential via (W(C)) ~ 
e ~ ltv QQ( L ) . This is what we call the "heavy quark potential". 

The expectation value of W(C) can be calculated at strong t'Hooft coupling in non- 
Abelian plasmas at large iV c that admit a weakly coupled dual gravity description according 
to AdS/CFT fl. More specifically, 

(W(C)) C FT = Z strmg (3) 

where Z string is the full supersymmetric string generating functional, which is defined in a 
10 dimensional background spacetime and includes a sum over all the string worldsheets 
whose boundary coincide with C. In the supergravity approximation A = g 2 N c ^> 1 and 
N c —>■ oo and, in this case, an infinitely massive excitation in the fundamental representation 
of SU(N C ) in the CFT is dual to a classical string in the bulk hanging down from a probe 
brane at infinity j^, 5]. Within this approximation Z str i n g ~ e* 5jvG and the dynamics of the 
string is given by the classical Nambu-Goto (NG) action (we neglect the contribution from 
other background fields such as the dilaton) 



S 



NG 



^ J d 2 a V-det h ab (4) 
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where h a b = G ^daX^dbX" (a, b = 1, 2), G^ u is the background bulk metric, a a = (r, a) are 
the internal world sheet coordinates, and = X^ir, a) is the embedding of the string in 



the 10- dimensional spacetime. For M = 4 SYM, the configuration that minimizes the action 
is a U-shaped curve that connects the string endpoints at the boundary and has a minimum 
at some U* in AdS 5 {4, 5]. 
The induced metric is 

deth ab = X' 2 -X 2 -(X-X') 2 (5) 

where X'^(t, a) = d a X^(r,a) and X M (r, cr) = d T X^(r, a). We choose a gauge where the 
coordinates of the static string are X^ = (t,x,0,0,U(x)), where r = t and a — x. We 
neglect the string dynamics in the 5 dimensional compact space and perform the calculation 
in real time at the boundary. In fact, fixing the extremal configuration in such a way implies 
the t — > 00 limit; while a Wick rotation is of course still possible (by switching to Euclidean 
metric, which would still give a complex expectation value for the Wilson loop) one can 
no longer perform an analytic continuation to imaginary time where the expectation value 
should be real. 
In this gauge, 

Sng = -7^1 [ dx^U' 2 + V{U) . (6) 
2na' J_ L/2 

where T — > 00 is the total time interval. Note that we have assumed that GqqGuu — 1 (which 
is in general valid when A — > 00) and defined V{U) = GoqG xx , which satisfies V(U) > for 
U G [Uh, 00). The equations of motion obtained from Eq. determine the classical string 
profile U c (x) as discussed in detail in refs. . The solution x = x(U c ) satisfies the 

following boundary condition 

-1/2 



which is used to obtain U» = U*(L, T). The expectation value of the Wilson loop is obtained 
by substituting the classical string solution U c (x) into the action in Eq. (jSJ). In general, when 
LT 1 the dominant contribution to the potential comes from the extremal worldsheet 
configuration described above and the potential is computed as a series in LT. Other 
configurations are expected to contribute significantly when LT > 1 [1JJ] . However, L < 1 /T 
is in fact the region of interest for bound states of very heavy quarks which have small radii. 

The heavy quark potential in the vacuum of TV = 4 SYM has the following simple 
analytical form (after subtracting the self-energy contribution from the infinitely massive 
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quarks) j^] 

Vqq( l ) = ~r( 1 / 4 )4X > ( 8 ) 
which may be compared to the standard SU(iV c ) Coulomb potential at large N c correspond- 
ing to weak coupling: 

Vcoui(L) = --L^ . (9) 
Sir L 

We now take into account thermal fluctuations around the classical solution U c (x). We 
shall show that the Wilson loop develops an imaginary part due to fluctuations near the 
bottom of the classical string configuration £/*. We consider long wavelength fluctuations of 
the string profile U c (x) — > U c (x)+5U(x) (with 5U' — > 0), which give the leading contribution 
to the string partition function in the supergravity approximation as follows 

Z string ~ J VX»e lSNG(x) ~ J V5U(x) e iS NG{u c +w) _ ^ 

When the bottom of the classical string is sufficiently close to the horizon (though still 
above it), the worldsheet fluctuations, SU(x), near x = (where U' c = 0) can change the 
overall sign of the argument of the NG square root and generate an imaginary contribution 
to the action. In this case, both U' c 2 and V(U C ) are small and the NG square root cannot 
be expanded in powers of 8U. 

The integral over SU(x) is performed by dividing the x interval in 2N parts such that 
xn = L/2, x_at = —L/2, Xj = jAx: 

Z strmg ~ J d(6U- N ) ■ ■ ■ d{5U N ) e'^ VW^W) . (ii) 

Near x = we expand U c {xj) ~ U* + x) U^(0)/2 and XJf + V{Uj) ~ C x x) + C 2 , with 

Ci = ^'(0) Uu'^O) + VI + 6U V;' + \bU 2 VI" 

* \u'M {2U'M + vl) > o 

C 2 = V* + 5UVl + Uu 2 V: . (12) 

Here, V* = V(U*), Vl = V'(U*) and so on. The imaginary part of the QQ potential arises 
from the region of 5U where C 2 < 0. 

We isolate the contribution to the path integral from x = Xj, 

Ij = / d{SU 3 ) exp -%— T {x)C x + C 2 ) 1/2 (13) 



where 5Uj m i n ^ max < are defined as the zeros of the argument of the square root. In this 
range of SUj the square root in the exponent develops an imaginary part. (The complement 
of the integration region in Eq. (1131) provides a correction to the real part of the potential 
due to fluctuations about the extremal configuration that is not considered here.) Note that 
in this case U* + 5U ~ Uh, i.e., the bottom of the fluctuating string touches the horizon. 
This may be viewed in the dual gauge theory as a process analogous to Landau damping 
of the static color fields which bind the quarks together, leading to the formation of two 
unbound heavy quarks in the high-temperature plasma. 

On the other hand, at lower temperatures on the order of the QCD crossover temperature 
and below, the dominant process should instead correspond to the breakup QQ — » (Qq) (Qq) 
into color-singlet heavy-light bound states; q stands for a light quark. Such tunnelling 
processes provide an exponentially small thermal width of deeply bound states 12]. It 



should be clear that the problem of quark tunneling cannot be solved rigorously since it 
involves genuinely non-perturbative QCD dynamics. However, the large mass of the heavy 



quark allows one to use the quasiclassical approximation [12|]. For a calculation of QQ — > 
(Qq) (Qq) via the AdS/CFT correspondence see Ref. Here, too, the temperature must 
be sufficiently low to allow for the formation of two new heavy-light quark bound states. In 
the gravity description, the final state corresponds to two strings connecting the Q-brane 
with the q-brane while in our approach they connect to the black-hole horizon. 

The factor 1/a! in the exponent implies that the leading order contribution is of order vA. 
In the supergravity approximation A ^> 1 and, thus, Ij can be computed in the saddle-point 
approximation. This gives 5U = —VJ./V" and so 



exp{-iTV QQ } = Y[lj 

j 

~exp(--^ 7 f/ dx v /-x*C 1 -V* + V: 2 /2V:' 

I 27r<y U\ x \ <Xc 

+ i j dx ^U' C 2 + V(U C ) \ (14) 

J\x\>Xc J J 



where x c = yJ—V* + VI 2 /2V")/Ci if this root is real, and x c = otherwise. The second 
contribution in the exponent gives the real part of the QQ potential which we drop from 
now on (see refs. 0, O,^, E|)- Performing the integral over \x\ < x c we find 

" vi v: 



Im V, 



QQ 



a' 



2K" Kj 



(15) 



where we used that U"(0) = VI/ 2. For J\f = 4 SYM at large A the potential entering the NG 
action is given by V(U) = (U 4 — U%) /R A , where R is the radius of AdSs and Uh = itR 2 T. 
This gives 

^V QQ = -^T*f±, (16) 

where a/A = R 2 /a' and ( = Uh/U* < 1. Note that the equation above applies only when 
C > ~ 0.76; otherwise Im Vqq = because the solution for x c from above ceases to 

exist. Thus, in the vicinity of this point the width generated by the imaginary part of the 
potential (see below) is small compared to the binding energy. 

The dependence of ( on L and T can be found from Eq. (17T) . At small LT we have 
LT = b( with b = 2r(3/4)/ v /^r(l/4) « 0.38. On the other hand, when the bottom of the 
classical string comes too close to the horizon, ( ~ 0.85, the U-shaped configuration used 



here receives higher-order corrections [ll| and cannot be used anymore. With ( ~ LT, the 
imaginary part of the potential (TlB"]) is smaller than the dominant contribution to the real 
part, eq. ([8]), by a factor ~ (LT) 4 . Here, we consider only temperatures such that thermal 
screening corrections to the real part of the potential are small; the bound state then probes 
the potential only in the region LT < 1. 

The imaginary part in Eq. ( TlBl) shifts the Bohr energy level obtained with the Coulomb- 
like vacuum potential ([8]), E — > E — iT; to first order, 



ny/X b 



48V2a 



45 



a T x 4 



(17) 



where denotes the unperturbed Coulomb ground state wave function and a = 
r(l/4) 4 /27r 2 \/AmQ is the Bohr radius. The width decreases with the quark mass and with 
the t' Hooft coupling, approximately as Tqq ~ l/Am^; it increases rapidly with the tem- 
perature, ~ T 4 . For m Q = 4.7 GeV, T = 0.3 GeV, V\ = 3 Q we obtain T T ~ 48 MeV. 

It is interesting to compare this result to the imaginary part of the heavy quark potential 
computed in pQCD: Im Vqq ~ — a 2 CFN c T(LT) 2 log (LT)' 1 [10]. This is smaller than 
R- e Vqq ~ a s Cp/L by one power of the coupling and three powers of LT. With L ~ 
(a s CVm<3) _1 the width decreases less rapidly with the quark mass than for f\f — 4 SYM at 
strong coupling. However, the numerical value of is on the order of tens of MeV, similar 
to what we obtain here. 
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We suggest that the width computed above is accessible experimentally through the 
suppression of T — > £ + £~ processes in heavy- ion collisions at RHIC or LHC. Neglecting 
"regeneration" of bound states from b and b quarks in the medium we can estimate the 
number of T mesons in the plasma at mid-rapidity, which have not decayed into unbound b 
and b quarks up to time t after the collision, from 

^ = -T T (T(t))N(t) 
->JV(t) ~ iV exp^- J dtV r {t) \ . (18) 

This solution assumes that Ty(T(t)) is a slowly varying function of time. The initial number 
of T states may be estimated from the multiplicity in p+p collisions times the number 
of binary collisions at a given impact parameter: N ~ N co uNj p . Thus, the integrated 
"nuclear modification factor" Raa for the process T — > £ + £~ is approximately given by 
RaaC^ ~^ £ + £~) — exp(— fx t), where F denotes a suitable average of T(T) over the lifetime 
of the quark-gluon plasma. Due to the strong temperature dependence of the width, this 
average is dominated by the early stage and thus we expect that F provides an estimate of 
the initial temperature in heavy- ion collisions via Eq. (TT7T) . For t = 5 fm/c and Ft = 48 MeV 
we obtain Raa — 0.3. An experimental estimate for the thermal quarkonium decay rate F 
could be obtained once a statistically significant detection of the T — > £ + £~ process has been 



achieved [151 ]. 
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